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We report the exact solution of the spectral problem for a graphene sheet framed by two armchair- 
and two zigzag-shaped boundaries. The solution is found for the ir electron Hamiltonian and gives, 
in particular, a closed analytic expression of edge-state energies in graphene. It is shown that the 
lower symmetry of graphene, in comparison with Cqh of 2D graphite, has a profound effect on the 
graphene band structure. This and other obtained results have far-reaching implications for the 
understanding of graphene electronics. Some of them are briefly discussed. 



Introduction. — Most theoretical studies on the elec- 
tronic properties ofgraphene start from either a Hiickel- 
type Hamiltonian EIJ3J1J3, S| , or k/kp versions of the 
Dirac Hamiltonian [1^311)111 > f° r a review, see [l(J. In 
one way or another, references are made to the symmetry 
points of 2D graphite band structure [ll| . At these points 
the valence and conduction 7r electron bands join each 
other, and the dispersion of electrons and holes is linear 
up to energies < 1 eV. This and the periodicity of two 
sublattices of 2D graphite unites electrons and holes near 
the Fermi energy with massless fermions. However, for 
finite-size graphene structures such as graphene ribbons, 
the structure illustrated in Fig. 1 is more relevant as a 
reference model than the 2D graphite lattice. Shown here 
is a plane macromolecule consisting of N X J\f hexagons 
which are arranged in sequences of N A/"-long oligomers 
of polyacene and coupled to each other via N C-C cova- 
lent bonds. All dangling bonds along the graphene edges 
are filled by hydrogen atoms. 

This model has been the focus of a number of works, 
but until now only approximate analytical solutions of 
the spectral problem have been proposed P, [H, @]. To 
begin, we present the exact description of the graphene 
7r electronic structure. It is shown that the spectrum is 
fully determined by the dispersion relation (which is dif- 
ferent from that known for 2D-graphite) supplemented 
by the generalized Lennard- Jones equation. This part 
of the description is substantially based on our previous 
studies of the band structure of conjugated oligomers 
flil . Ha |. Next, we discuss applications of the obtained 
equations in the context of graphene and its daughter lat- 
tices, armchair and zigzag graphene ribbons (GRs) and 
carbon nanotubes (CNTs). In conclusion, we express the 
edge-state spectrum in terms of elementary functions and 
show its agreement with the exact results. In what fol- 
lows, C-C hopping integral t is the only parameter of 
the semi-empirical Hamiltonian of the graphene macro- 
molecule. This parameter is used as a unit of energy. 

Exact solution of the eigenvalue problem for 
graphene. — By taking an appropriate representa- 
tion of the molecular orbitals of graphene, we reduce 
the initial two-dimensional Schrodinger problem to 
M independent sets of 27V one-dimensional equations. 
Each set describes a hypothetical oligomer consisting 
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FIG. 1: Indication of labels of carbon atoms used in the 
present description of the n electron spectrum of N x TV hon- 
eycomb lattice. Right inset shows symmetry points of 2D 
graphite (hexagon vertices) and graphene (crosses). In lower 
part, dashed- framed block (oligomer of polyacene with the 
length aj\f) can be thought as a monomer M in an M-oligomer. 
The monomer Green's function matrix elements completely 
determine electron spectra of M-oligomers [12l.ll3{. InEq. ([3|, 
the role of these matrix elements is played by gj r and g\ v 



of N monomers, as is illustrated in Fig. 1 by dashed, 
numbered frames. 

Exploiting m, n, and a — l,X,p,r labeling explained 
in Fig. 1, the ir electron wave function that satisfies the 
Schrodinger equation H 1 ^ — with the tight-binding 
Hamiltonian of Nx J\f graphene sheet, can be represented 
as follows 

* = EE E ^ 

m,n,a \m,n,a), (1) 

m n a=l,r,X,p 

where |m, n, a) is the 2p z orbital at the ath atom of ben- 
zene ring with coordinates {m, n} with a summation run- 
ning over all sites of the honeycomb lattice (|m, n, a) = 
if these lattice sites are empty, for example, 1 1, N+l, a) — 
0), and 
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A/", and coefficients 0£ a , 
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In the latter expansion, a? n 
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#,j> 9i,i = 9lr, V jff > r = 4cos^(0/2), 
E[E 2 - 1 - 4cos 2 (^/2)], and zeros of Vj = 
[E 2 — 4cos 2 (£j/2)] 2 — E 2 determine the spectrum of an 
N- long acene. 

The set of equations is central in this otherwise 
standard derivation. As already mentioned, it appears in 
the theory of M-oligomers, M-M-. . . -M, where the energy 
dependent quantities of gj l and ; are the monomer 
Green's function matrix elements referring to the same 
(left or right) or different binding atoms of a monomer M; 
see Fig. 1. This analogy was first noticed by Klymenko 

Finding the eigenvalues of ([3]) and thus solving the 
eigenvalue problem dct(H — EI) = yields 



cosk = f{gli,gfi r ) 



1 



Mr 



1 



(a;.,) 2 (ah ' 



(4) 



where n and E are interrelated via 
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Formally the same equations as (j4|) and (l5l) appea r in the 
tight-binding description of M-oligomers [TJ,ll3|. As a 
particular case, in the Lennard- Jones theory of polyenes 
(M = C=C) the right hand side of Eq. © is inde- 
pendent of energy and equal to a certain constant. 

By taking into account the explicit expressions of gj t 
and gi l , Eqs. (H|) and can be transformed into 
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respectively. Within the interval < k < tt, Eq. ([7]) 
with sign plus or minus has N solutions which determine 
j-dependent quantum numbers Kj „, v = 0,1, . ■ ■ , N — 1 
and hence, the spectrum of the graphene sheet, E = 
±E\ .. To be precise, Eqs. © and {7]) determine ANN 

of the total number 2N(2N + 1) of the tt electron levels. 
Additionally, there are two A^-fold degenerate levels with 
energies ±1. These can be proven to be the states with 
zero wave- function amplitudes at the I and r sites, thus 
making them of no interest here. Below, only ANN tt 
electron states, \j, v), are considered. 

Spectra of graphene daughter lattices. — Equation (|6|) 
remains valid for periodic boundary conditions (PBCs), 



in which case the tt electron spectrum is fully deter- 
mined by this single equation, where n — k; — 2ttI/N, 
I = 0,1,2,. . . ,N-1 and ^ = 2nj/N, j = 0,1,2,. . . J\f-1. 
Thus defined, the dispersion relation reads 



[^(fer, ky)] 2 = 1 ± 4 cos(ak y /2) cos(V3ak x /2) 
+ A cos 2 (ak y / '2), 



(8) 



where the minimal translation distance a is indicated in 
Fig. 1; the correspondence between the continuos vari- 
ables in the dispersion relation and discrete quantum 
numbers in Eq. ^ is as follows: ^/3ak x <-*■ k, ak y <-» £j. 

For the system in focus, the usage of Eq. {8]| should be 
restricted to the range < v3|^k|) \k y \ < ir/a. It is easy 
to see that, within this range, there are only two points, 
k x = Q,ak y = 2n/3 and k x — 0,ak y = —2tt/3, where 
E^(k x , k y ) = (instead of six for 2D graphite |ll|). Not 
far away from these points, the dispersion relation ([8]) 
can be approximated by the familiar formula 



(k x , k v ) = ±^sjkl + {\ky\ - 2-K/Zaf 



(9) 



that is by a linear form of dispersion, if it is expressed 
in terms of the deviation of k from the points of zero 
energy. The reduced number of these zero points is a 
direct consequence of the lower symmetry in graphene in 
comparison with 2D graphite. 

Equations (|6|) and ([7|) can be easily used for an instruc- 
tive description of graphene daughter lattices, specifi- 
cally, armchair and zigzag GRs and CNTs. However, the 
cases of finite and infinite (i.e., independent of the bound- 
ary conditions) systems must be clearly distinguished. 
For infinite armchair GRs, the spectrum is completely 
determined by Eq. where —tt < n — \fiak x < tt. For 
infinite zigzag GRs, it is £ = ak y that should be treated 



E ±z = 1±4|cos(&/2)cos(k/2)| + 4cos 2 (^/2), (6) q 



as a continuous variable in the two equations, 



and 



Formal description of armchair and zigzag CNT spec- 
tra is exactly the same, as that of zigzag and armchair 
GRs, respectively. The only difference is that the dis- 
crete quantum number is not determined by the open 
boundary conditions [as in Eq. |2])] but by the PBCs. 
As a result, the use of the dispersion relation ([8]) ex- 
tends to the range < v^l k x \ < Tr/a, < \ky\ 5; 2tt/cl 
and < \/3|fcx| < 2n/a, < \k y \ < n/a, for zigzag 
and armchair CNTs, respectively. This means that the 
band structure of zigzag CNTs, as compared with arm- 
chair GRs, has two new points of zero energy, E^(k x = 
0,k y — ±An/3a) = 0. In armchair CNTs, four special 
points appear, (k x = ±27r/\/3a, k y — ±2n/3a). Further 
discussion of CNTs and GRs spectra can be found else- 
where (ltij . 

Some of the essential parts of the above discussion are 
exemplified in Fig. 2. Represented in its upper part is 
the tt electron spectrum of a 21 x 21 graphene sheet, where 



FIG. 2: Graphene spectrum E + , E > 0, according Eqs. 1SJ and (J7J). Squares and circles correspond to signs plus and minus 

in these equations, as explained in the text; different colors indicate different values of j. From left to right: General view, 
projections £ = 0, and k — 0. Af = N = 21 and A/" = 75, N = 60 for upper and lower panels, respectively. Imaginary quantum 
values of k are shown on 5 continuation of k axis. The corresponding energy levels (edge-state levels) are seen as dots (upper 
panels) and black line (lower panels) in 5-£ plane. 



each value of j has its own color and circles and squares 
correspond to "minus" and "plus" branches, respectively, 
of Eq. ([6]). Three panels show (from left to right) energies 
of j, v levels, and crossections of the spectrum by planes 
k-E and £-E. Levels which are represented by circles 
in 5-£ plane correspond to imaginary values of k. They 
associate with electron states which are localized near 
zigzag-shaped boundaries. The conjugated part of the 
spectrum, j, v levels with negative energies, is the mirror 
reflection in k-£ and planes. 

Transformation of the spectrum with the increase of N 
and Af is illustrated by the lower panels in Fig. 2. These 
give a visual representation of the graphene band struc- 
ture E ± (k x , k y ). Note that the spectra shown in the mid 
and right lower panels have the same appearance as those 
which have been obtained in computational modeling of 
armchair and zigzag GRs [l[. 

Spectrum of graphene edge states. — By denoting that 
£j — 2ir/3 = qj , it can be shown that Eq. ([7]) has N real 
solutions, nj < kJx<. . . <nJ N _ 2 < n~ N _ l7 if qj < q c , 
where q c can be found from Eq qC = ±[1 — 2cos(7r/3 + 
q c /2)\. If TV » 1 then, q c = (vtiV)" 1 . The smallest of 
the solutions for k becomes imaginary if qj > q c : kJ = 
iSj, j > J*, where j* = [2(7V + l)/3]; [A] denotes a 
minimal integer of rational number A. Note that kJ* q = 
Sj* = 0, if (Af + l)/3 is an integer. In Fig. 2, imaginary 



values of k„^„„ n = iSj are shown on an extension of the 

are very close but are 



k axis. Energy levels E~~_ 



never equal to zero; see below. 

The energies which satisfy Eqs. © and 0, and fall 
into the interval -(27V + l)" 1 < E < (27V + l) -1 , cor- 
respond to imaginary values of kJ and hence, to elec- 
tron states, decaying towards the mid of the graphene 
sheet along the armchair direction. These states (a kind 
of Tamm/Shockley surface states in molecular structures 
[l7j |) have been discussed by many authors in the con- 
text of zigzag graphene ribbons [H, H, 0, H, 0] j where such 
states associate with decaying modes of the transverse 
electron motion. In distinction from the Tamm states, 
edge states of graphene are not fully localized: They de- 
cay in the (transverse) armchair direction. In the (lon- 
gitudinal) zigzag direction, these states are delocalized 
and can be described by a superposition of propagating 
states. 

By restricting ourselves to imaginary values of kJ — 
iSj, that is by the spectrum of decaying modes, we can 
rewrite the minus branch of Eq. (JSJ) in the form 
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FIG. 3: Energy of edge states, E^.,, as a function of 
Exact dependence and Eq. fjl2p are represented by a single 
solid line; approximation (|13|l for SjN » 1 and n/3 — qj « 
1 are shown by dashed and dotted lines, respectively. 



the Fermi energy. This result devalues the concept of a 
"zero mode" with no dispersion. Secondarily, this solves 
a number of long standing problems, which have been 
the subject of a considerable computational and analyti- 
cal effort. Not immediate but a straightforward applica- 
tion of the obtained solution is an accurate description 
of the spectrum of achiral graphene ribbons and carbon 
nanotubes near the point of neutrality [l6[. Altogether 
this forms a new platform for arguable interpretation and 
modeling of the electronic properties of graphene and its 
daughter structures. 
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where Sj must satisfy 
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The two equations above give the exact position of 
edge-state levels in the graphene spectrum. An approxi- 
mate solution of Eqs. (fTTj) , Sj = —2 In [2 cos(£j/2)], yields 
the values of edge-state energies, 
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(12) 

which even in the logarithmic scale are indistinguishable 
from the exact solution; see Fig. 3. For values of qj, 
which are close to 7r/3 and under the condition SjN » 
1, Eq. (HU simplifies to 
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j>j*= Sjexp(-SjN), 



S 3 N » 1, 



Tr-nj/(M+l)] 2N ,7r/3-<& « 1.(13) 



The exponential behavior shown here has been ob- 
served previously in many numerical models, e.g., [l|, 0, 
5], and an analytical description was given by Brey and 
Fertig Q. Qualitatively, the latter agrees with Eq. (fTB")) . 
but the functional form of the exponential factor is very 
different from our exact result. Further comments on this 
point can be found in Ref. [181 ]. 

To summarize our findings, we have presented an ex- 
act quantitative description of a rectangular sheet of 
graphene which brings to light the crucial role of zigzag 
boundaries in determining the tt electron spectrum near 
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